Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
竹田 武司
JAEA-Data/Code 2023-012, 75 Pages, 2023/10
JAEA-Data-Code-2023-012.pdf:4.45MB
ROSA-V計画において、大型非定常実験装置(LSTF)を用いた実験(実験番号: TR-LF-15)が2014年6月11日に行われた。ROSA/LSTFTR-LF-15実験では、加圧水型原子炉(PWR)のポンプシール冷却材喪失事故(LOCA)を伴う、補助給水機能喪失を特徴とするTMLB'のシナリオでの全交流電源喪失時のアクシデントマネジメント(AM)策を模擬した。ポンプシールLOCAは、0.1%低温側配管破断により模擬した。このとき、非常用炉心冷却系(ECCS)である高圧注入系及び低圧注入系の全故障とともに、ECCSの蓄圧注入タンクから一次系への非凝縮性ガス(窒素ガス)の流入を仮定した。蒸気発生器(SG)二次側水位が特定の低水位まで低下すると、一次系圧力は上昇に転じた。SG二次側水位喪失後、加圧器の安全弁が周期的に開いたため、一次冷却材の喪失につながった。故に、高圧条件でボイルオフによる炉心露出が生じた。模擬燃料棒被覆管表面温度の10Kの上昇を確認した時点で、SG二次側減圧を一番目のAM策として開始した。このAM策では、両SGの安全弁を開放した。また、一番目のAM策開始後少し遅れた時点で、加圧器の安全弁の開放による一次系減圧を二番目のAM策として開始した。さらに、一番目のAM策に従いSG二次側圧力が1.0MPaに低下した時点で、低水頭ポンプによる給水ラインから両SG二次側への注水を三番目のAM策として開始した。三番目のAM策の開始直後、SG二次系からの除熱が再開したため、一次系圧力の低下が促進された。蓄圧注入系から両低温側配管への冷却材注入による炉心水位の回復により、全炉心はクエンチした。窒素ガスがSGU字管内に蓄積したため、一次系の減圧率は低下した。本報告書は、ROSA/LSTFTR-LF-15実験の手順、条件および実験で観察された主な結果をまとめたものである。
中村 いずみ*; 大谷 章仁*; 奥田 幸彦; 渡壁 智祥; 滝藤 聖崇; 奥田 貴大; 嶋津 龍弥*; 酒井 理哉*; 渋谷 忠弘*; 白鳥 正樹*
第10回構造物の安全性・信頼性に関する国内シンポジウム(JCOSSAR2023)講演論文集(インターネット), p.143 - 149, 2023/10
原子力発電施設における配管系の耐震設計では、設計対象を弾性はり要素でモデル化し、弾性解析に基づく保守的な応力評価を実施している。一方、これまでに実施された多数の実験結果から、配管系は設計の想定を超えるような地震入力下では弾塑性挙動を示し、破損に至るまでには大きな裕度を有していると認識されている。このような状況を踏まえ、適切な保守性と合理性を有する耐震評価のため、弾塑性応答挙動を考慮した新たな耐震設計・評価手法の構築を目指し、2019年に日本機械学会より発電用原子力設備規格設計・建設規格の事例規格が発刊された。初版発刊後は事例規格の継続的な改善のために議論と検討を進め、2022年には疲労評価に用いるサイクルカウント法等に修正を加えた改訂版の発刊が決定した。また、次期改訂に向け、配管支持構造物の弾塑性評価を規格に取り入れる議論が進められている。本稿では、2022年の事例規格における主要な改訂内容、改訂の背景、次期改訂に向けた取り組み状況及び今後の課題について紹介する。
Zhang, H.*; 梅原 裕太郎*; 吉田 啓之; 森 昌司*
International Journal of Heat and Mass Transfer, 211, p.124253_1 - 124253_13, 2023/09
被引用回数：2 パーセンタイル：59.37(Thermodynamics)In this study, the characteristics of disturbance waves, including velocity, longitudinal size, and frequency, are comprehensively investigated by conducting gas-liquid upward annular flow experiments in a 5 mm tube through the constant electric current method (CECM). The film thickness time trace is recorded concerning time, and the information on disturbance waves is subsequently extracted. The effect of liquid and gas flow rate, density ratio, and surface tension on the disturbance wave is quantitatively analyzed. A predictive model of the wave velocity based on the shear stress balance of the liquid film is proposed. A physical model for wave frequency prediction is derived on the basis of the mass conservation equation. A simple empirical model with good prediction accuracy of wave frequency is also derived. Compared with existing correlations available in the literature, the newly derived models show better performance under a wide range of flow conditions.
有吉 玄; 猿田 晃一; 粉川 広行; 二川 正敏; 前野 航希*; Li, Y.*; 筒井 喜平*
Proceedings of 20th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-20) (Internet), p.1407 - 1420, 2023/08
水銀を核破砕標的とするパルス中性子源では、陽子線励起圧力波による標的容器のキャビテーション損傷が問題視されている。このような損傷の抑制手法として、旋回流式気泡注入器で水銀中に微小He気泡を混合し、その圧縮性効果を用いる手法や、標的容器内側に流路壁を追加した二重壁流路構造を設ける手法が採用されている。本研究では、二重壁流路体系において内壁部に生じ得るキャビテーション由来の貫通損傷を標的容器外部から検出する異常診断技術の確立を目指す。本論文では、貫通損傷が生じた場合の水銀流れを有限要素法に基づく流動解析で明らかにし、流体励起・音響振動の観点から、損傷の影響を評価した結果を報告する。
澤口 拓磨; 高井 静霞; 笹川 剛; 打越 絵美子*; 嶋 洋佑*; 武田 聖司
MRS Advances (Internet), 8(6), p.243 - 249, 2023/06
放射性廃棄物の中深度処分では、モニタリング用のボーリング孔内が適切に閉塞されているかを確認するための手法を事前に整備しておく必要がある。そこで本研究では、堆積岩地域を想定し、どのような埋戻し設計条件であればボーリング孔内が有意な移行経路とならないかを把握し、ボーリング孔閉塞に係る確認ポイントを明らかにするため、埋戻されたボーリング孔を有する水理地質構造に対する地下水流動解析を実施した。その結果、ボーリング孔や掘削損傷領域(BDZ)が移行経路とならないための条件として、ベントナイト系材料の透水係数を母岩と同等以下にすること、BDZにグラウトを充填することなどが示された。
Zhang, H.*; 森 昌司*; 久野 努*; 吉田 啓之
International Journal of Multiphase Flow, 159, p.104342_1 - 104342_15, 2023/02
被引用回数：5 パーセンタイル：65.31(Mechanics)At the gas-liquid interface, disturbance waves with various heights commonly arise. Because the thin liquid film between two successive disturbance waves leads to the dryout on the heating surface and limits the performance of the heating components, complete knowledge of the film thickness and disturbance waves is crucial. However, most existing experimental data and analysis are limited to near atmospheric conditions, and the liquid gas density ratio is much higher than in practical nuclear applications. In this study, we use nitrogen gas and water under the system pressures of 0.2 and 0.4 MPa and HFC134a gas and water under the system pressure of 0.7 MPa. The density ratio can be varied between 32 and 434. Based on experimental results, two models are proposed to predict the averaged film thickness and wave height, which are then compared with previous data.
渡辺 勢也*; 河原 淳*; 青木 尊之*; 杉原 健太; 高瀬 慎介*; 森口 周二*; 橋本 博公*
Engineering Applications of Computational Fluid Mechanics, 17(1), p.2211143_1 - 2211143_23, 2023/00
被引用回数：1 パーセンタイル：61.52(Engineering, Multidisciplinary)津波氾濫や大雨の斜面災害では、多くの浮遊物や流木、丸太などが流れに含まれる。剛体衝突による構造物の被害は、水圧による被害よりも深刻である。浮遊物体を含む自由表面流れを研究するためには、大規模計算が可能で高性能な自由表面流のシミュレーションコードの開発が必要となる。本研究ではキュムラント格子ボルツマン法と粒子ベースの剛体シミュレーションを組み合わせた単相自由表面モデルを提案する。剛体間の接触相互作用は個別要素法で計算される。解析精度の向上と計算の高速化のために8分木ベースの局所細分化格子法を導入し、自由表面や固体表面近傍には高解像度の格子を割り当てた。提案モデルの精度検証のために、八戸工業大学と神戸大学の15m水槽と70m水槽で2種類の津波実験を行った。シミュレーションの結果、漂流速度、捕捉木片数、積層角について実験と良い一致を示した。
吉田 啓之; 堀口 直樹; 小野 綾子; 古市 肇*; 上遠野 健一*
Proceedings of 29th International Conference on Nuclear Engineering (ICONE 29) (Internet), 7 Pages, 2022/08
About the boiling transition (BT) that determines the maximum thermal output of the BWR, it is considered that the spacers have significant effects on the occurrence of BT. And occurrence conditions of BT can be changed by devising the spacer shapes. In the light water cooled fast reactor: RBWR, thermal-hydraulics conditions are more severe than the current BWR. Then, the effect of the spacer on BT should be sufficiently utilized in the RBWR. In the thermal-hydraulics design for the current BWR, large-scale tests were carried out and used to evaluate BT conditions. The RBWR is still in the design stage, and there is room to be changed to many parameters. Then, it is not reasonable to determine the shape of the spacer by evaluation only for large-scale tests. On the other hand, by applying a two-phase CFD method with remarkable development in recent years, we can develop a model that can predict the effect of spacers mechanistically. This research used the detailed two-phase flow simulation code TPFIT developed by JAEA to simulate annular dispersed flow in RBWR subchannels. In the occurrence of BT, it is considered that the two-phase flow pattern is the annular dispersed flow, and we want to evaluate the effects of spacer shape on annular dispersed flow in RBWR subchannels. As the first step of this research, we performed numerical simulations of annular dispersed flow in the simplified subchannel of RBWR. We used a circular tube with the same hydraulic diameter as the RBWR subchannel to consider the basic effects of spacer on the annular dispersed flow. As a simulation parameter, we choose the existence of the spacer. The spacer used in the simulation has a simplified shape and the same blockage ratio as the RBWR. In this paper, we describe the result of numerical simulation. We evaluated droplets' size and velocity based on simulation results for the spacer's existence and non-existence cases.
Bateman, K.*; 村山 翔太*; 花町 優次*; Wilson, J.*; 瀬田 孝将*; 天野 由記; 久保田 満*; 大内 祐司*; 舘 幸男
Minerals (Internet), 12(7), p.883_1 - 883_20, 2022/07
被引用回数：0 パーセンタイル：0.02(Geochemistry & Geophysics)The construction of a repository for geological disposal of radioactive waste will include the use of cement-based materials. Following closure, groundwater will saturate the repository and the extensive use of cement will result in the development of a highly alkaline porewater, pH 
12.5; this fluid will migrate into and react with the host rock. The chemistry of the fluid will evolve over time, initially high [Na] and [K], evolving to a Ca-rich fluid, and finally returning to the groundwater composition. This evolving chemistry will affect the long-term performance of the repository, altering the physical and chemical properties, including radionuclide behaviour. Understanding these changes forms the basis for predicting the long-term evolution of the repository. This study focused on the determination of the nature and extent of the chemical reaction, as well as the formation and persistence of secondary mineral phases within a granite, comparing data from sequential flow experiments with the results of reactive transport modelling. The reaction of the granite with the cement leachates resulted in small changes in pH and the precipitation of calcium aluminum silicate hydrate (C-(A-)S-H) phases of varying compositions, of greatest abundance with the Ca-rich fluid. As the system evolved, secondary C-(A-)S-H phases re-dissolved, partly replaced by zeolites. This general sequence was successfully simulated using reactive transport modelling.
鳥川 智旦*; 大平 直也*; 伊藤 大介*; 伊藤 啓*; 齊藤 泰司*; 松下 健太郎; 江連 俊樹; 田中 正暁
混相流, 36(1), p.63 - 69, 2022/03
ナトリウム冷却高速炉(SFR)において、炉容器上部プレナム内の自由界面で、自由表面渦によるカバーガス巻込みが発生し、巻き込まれたガスが炉心を通過することで炉出力の擾乱を生じる可能性がある。そのため、巻き込みガスの流量を正確に予測するための評価モデルの開発が重要となる。既往研究における単一渦のガス巻込み試験では、一般的に、ガス巻込みは上部タンクの自由表面渦によって引き起こされ、吸込み管によって下部タンクにおいてガスが液体から分離される。しかし、これらの研究では、上部タンクと下部タンクとの間の圧力差の影響に着目していない。本研究では、上下部タンク間の圧力差の影響に着目した単一渦のガス巻込み試験を行った。上部タンクと下部タンクとの圧力差は、下部タンクのガス圧力を変更して制御した。その結果、上部タンクと下部タンクの圧力差が増加するにつれ、巻込みガス流量も増加することが分かった。また、吸込み管内の旋回環状流の可視化により、巻込みガス流量が増加すると吸込み管内の圧力降下が大きくなることが分かり、旋回環状流領域における圧力降下に基づいた評価モデルによってガス巻込み流量を予測できることが示唆される。
小野 綾子; 山下 晋; 鈴木 貴行*; 吉田 啓之
Proceedings of 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19) (Internet), 16 Pages, 2022/03
原子力機構では、軽水炉燃料の安全性評価において必須である限界熱流束の評価において、新型燃料設計にかかわる大型モックアップ試験によるコストの削減や、想定外事象に対応するためにモックアップ試験の試験範囲よりも幅広い範囲において、限界熱流束を機構論に基づき評価する研究に着手している。本研究では、機構論的流動解析手法であるJUPITERを用いて、スワールベーンおよびスプリットベーン付きの燃料集合体内の国際ベンチマーク問題を対象とした単相流解析および同体系における二相流解析を実施し、各種ベーンによる流動場および気泡挙動に与える効果、解析における課題の抽出を行った結果を報告する。
羽賀 勝洋; 粉川 広行; 直江 崇; 涌井 隆; 若井 栄一; 二川 正敏
Proceedings of 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19) (Internet), 13 Pages, 2022/03
J-PARCではMWクラスの核破砕中性子源を実現するために水銀ターゲットの流路構造としてクロスフロー型ターゲットを開発し、流路構造の改良を継続してきた。高出力の短パルス陽子ビームが水銀ターゲットに入射すると、急激な水銀の発熱と体積膨張により最大40MPaにも達する圧力波が誘起され、ターゲット容器にキャビテーション損傷を生じさせる。このため、水銀流路にマイクロバブル生成器を配置し、バブルの収縮により水銀の体積膨張をクッションのように吸収することで圧力波を低減する技術や、陽子ビームが入射する容器壁に内壁を設けて、狭隘流路に形成される速い水銀流れの大きな速度勾配を利用してキャビテーション損傷を低減する技術などを水銀ターゲット容器の流路構造に導入した。これらの技術開発により、2020年には36.5時間の1MW連続運転を成功させ、2021年4月から最大740kWの高出力で長期の安定な利用運転を達成した。本報告は、主に水銀ターゲット容器の熱流動設計に関して1MW運転を実現するまでの技術開発をまとめたものである。
Park, Y.-J.*; 澤田 淳; 小堤 健紀*; 田中 達也*; 橋本 秀爾*; 森田 豊*
Proceedings of 3rd International Conference on Discrete Fracture Network Engineering (DFNE 2022) (Internet), 8 Pages, 2022/00
高レベル放射性廃棄物の地層処分の安全評価には地層中の長期にわたる地下水流動と核種移行プロセスの把握が求められる。沿岸部地下環境において、地下水流動は海水起源の塩水と陸水起源の淡水の密度差による複雑な相互作用の影響を受ける。加えて、数百万年の長期においては、海進・海退に伴う海水準変動の影響を受ける。本研究では、そのような沿岸域における亀裂性の結晶質岩を対象とした処分場の地下水流動と核種移行を評価するため、塩分濃度と地下水流速などの地下水環境の長期的な変遷を評価するための広域スケールとブロックスケールを組み合わせた評価フレームを構築した。
大川 富雄*; 森 昌司*; Liu, W.*; 小瀬 裕男*; 吉田 啓之; 小野 綾子
日本原子力学会誌ATOMO
, 63(12), p.820 - 824, 2021/12
原子炉設計における効率的な燃料設計および最適な安全評価のために、機構論に基づいた限界熱流束評価技術が望まれている。この長年の技術課題は、近年の詳細解析技術及び計測技術の体系的統合を段階的に進めることで、打開できる可能性がある。このため本研究専門委員会では、将来的な限界熱流束評価技術の構築に向けて、過去の膨大な研究を精査することで必要な知見を整理する。これらの議論を通して、原子炉における機構論的限界熱流束評価技術に必要な研究課題を提示する。Part2では、これまでの限界熱流束機構に関する基礎研究や限界熱流束の予測手法確立に必須な数値解析手法の発展にふれ、課題提起を行う。
Han, X.*; Shen, X.*; 山本 俊弘*; 中島 健*; 孫 昊旻; 日引 俊*
International Journal of Heat and Mass Transfer, 178, p.121637_1 - 121637_24, 2021/10
被引用回数：7 パーセンタイル：58.99(Thermodynamics)This paper studies the flow regimes, their transitions and the drift-flux correlations in upward gas-liquid two-phase flows in vertical rod bundle flow channels. The flows are classified into 5 flow regimes, namely, bubbly, finely dispersed bubbly, cap-bubbly, churn and annular flows according to their different flow characteristics. Transition criteria between the flow regimes are proposed mechanistically. Those criteria can correctly predict 83% of the existing experimental observation of the flow regime. The drift-flux correlations for the distribution parameter and the drift velocity are also improved. The void fractions predicted by those correlations are compared with the existing experimental data, showing satisfactory agreement with mean relative error of 8%.
Bateman, K.; 村山 翔太*; 花町 優次*; Wilson, J.*; 瀬田 孝将*; 天野 由記; 久保田 満*; 大内 祐司*; 舘 幸男
Minerals (Internet), 11(9), p.1026_1 - 1026_23, 2021/09
被引用回数：2 パーセンタイル：23.36(Geochemistry & Geophysics)The construction of a repository for geological disposal of radioactive waste will include the use of cement-based materials. Following closure, groundwater will saturate the repository and the extensive use of cement will result in the development of a highly alkaline porewater, pH 12.5. This fluid will migrate into and react with the host rock. The chemistry of the fluid will evolve over time, initially high [Na] and [K], evolving to a Ca-rich fluid and finally returning to the groundwater composition. This evolving chemistry will affect the long-term performance of the repository altering the physical and chemical properties, including radionuclide behaviour. Understanding these changes forms the basis for predicting the long-term evolution of the repository. This study focused on the determination of the nature and extent of the chemical reaction; the formation and persistence of secondary mineral phases within an argillaceous mudstone, comparing both data from sequential flow experiments with the results of reactive transport modeling. The reaction of the mudstone with the cement leachates resulted in small changes in pH but the precipitation of calcium aluminium silicate hydrate (C-A-S-H) phases of varying compositions. With the change to the groundwater secondary C-(A-)S-H phases re-dissolved being replaced by secondary carbonates. This general sequence was successfully simulated by the reactive transport model simulations.
長尾 理那; 滑川 麻紀*; 戸塚 真義*; 仲田 久和; 坂井 章浩
JAEA-Technology 2021-009, 139 Pages, 2021/06
JAEA-Technology-2021-009.pdf:13.96MB
日本原子力研究開発機構は、低レベル放射性廃棄物のうち研究施設等廃棄物の埋設処分事業の実施主体となっている。低レベル放射性廃棄物を処分する方法としてコンクリートピット処分を検討している。コンクリートピット施設は、地下水位より深い場所に設計するため、地下水の流れとともに放射性核種が移行すると可能性があると考えられている。そのため、コンクリートピット処分施設の安全性を説明するために、地下水の流れやコンクリートピット施設からの浸出水量を調べる必要がある。そこで、本報告書では、有限要素法による二次元地下水流動解析コード(MIG2DF)を用いて、充填覆土の透水係数やベントナイト混合土の設計等の施設の設計条件を変えた感度解析を実施し、コンクリートピット施設からの浸出水量を算出した。また、長期にわたる管理期間中に発生の可能性があるベントナイト混合土の劣化について考慮した評価も行うこととした。解析の結果、ベントナイト混合土が健全な場合、側部覆土の透水係数を低くすることにより浸出水量が減少していた。このことから、側部覆土の透水係数を低くすることによる浸出水量の低減が可能であることが示唆された。しかし、コンクリートピット施設の上部のベントナイト混合土に劣化が生じた場合、側部覆土の透水性を低くしても、浸出水量の大幅な低減は見られなかった。そのため、コンクリートピット施設の上部のベントナイト混合土に劣化が生じる可能性を考慮し、側部覆土のみではなく、コンクリートピット施設の上部の覆土についても透水性の低い覆土を充填することを検討する必要があると考えられる。
Bateman, K.; 天野 由記; 久保田 満*; 大内 祐司*; 舘 幸男
Minerals (Internet), 11(6), p.588_1 - 588_19, 2021/06
被引用回数：4 パーセンタイル：45.1(Geochemistry & Geophysics)The alkaline fluid arising from the cement-based materials will react with host rock to form a chemically disturbed zone around the geological disposal repository. To understand how these conditions may evolve, a series of batch and flow experiments were conducted with Horonobe mudstone and fluids representative of the alkaline leachates expected from both Ordinary Portland Cement (OPC) and Low Alkali Cement (LAC) leachates. The impact of the LAC was more limited, compared to the OPC leachates. Ion exchange reactions had an impact initially but thereafter the reaction was dominated by primary mineral dissolution, and in the case of OPC leachates precipitation of secondary C-S-H phases. The flow experiments revealed that precipitation of the secondary phases was restricted to close to the initial contact zone of the fluids and mudstone. The experimental results demonstrate that a combination of both batch and flow-through experiments can provide the insights required for the understanding of the key geochemical interactions and the impact of transport, allowing the spatial as well as temporal evolution of the alkaline leachate / mudstone system to be determined.
竹田 武司
JAEA-Data/Code 2021-006, 61 Pages, 2021/04
JAEA-Data-Code-2021-006.pdf:2.78MB
ROSA-V計画において、大型非定常実験装置(LSTF)を用いた実験(実験番号: SB-PV-09)が2005年11月17日に行われた。ROSA/LSTF SB-PV-09実験では、加圧水型原子炉(PWR)の1.9%圧力容器頂部小破断冷却材喪失事故を模擬した。このとき、非常用炉心冷却系(ECCS)である高圧注入系の全故障と蓄圧注入(ACC)タンクから一次系への非凝縮性ガス(窒素ガス)の流入を仮定した。実験では、上部ヘッドに形成される水位が破断流量に影響を与えることを見出した。アクシデントマネジメント(AM)策として、両ループの蒸気発生器(SG)逃し弁開放によるSG二次側減圧を炉心出口最高温度が623Kに到達した時点で開始した。SG二次側圧力が一次系圧力に低下するまで、このAM策は一次系減圧に対して有効とならなかった。一方、炉心出口温度の応答が遅くかつ緩慢であるため、模擬燃料棒の被覆管表面最高温度がLSTFの炉心保護のために予め決定した値(958K)を超えたとき、炉心出力は自動的に低下した。炉心出力の自動低下後、低温側配管内でのACC水と蒸気の凝縮により両ループのループシールクリアリング(LSC)が誘発された。LSC後、炉心水位が回復して炉心はクエンチした。ACCタンクから窒素ガスの流入開始後、一次系とSG二次側の圧力差が大きくなった。ECCSである低圧注入系の作動を通じた継続的な炉心冷却を確認後、実験を終了した。本報告書は、ROSA/LSTF SB-PV-09実験の手順、条件および実験で観察された主な結果をまとめたものである。
岡垣 百合亜; 与能本 泰介; 石垣 将宏; 廣瀬 意育
Fluids (Internet), 6(2), p.80_1 - 80_17, 2021/02
Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in the convection term of the volume fraction equation. Thus, in this study, we focused on an interface compression (IC) method for such a VOF approach; this technique prevents numerical diffusion issues and maintains boundedness and conservation with negative diffusion. First, on a sufficiently high mesh resolution and without the IC method, the validation process was considered by comparing the amplitude growth of the interfacial wave between a two-dimensional gas sheet and a quiescent liquid using the linear theory. The disturbance growth rates were consistent with the linear theory, and the validation process was considered appropriate. Then, this validation process confirmed the effects of the IC method on numerical diffusion, and we derived the optimum value of the IC coefficient, which is the parameter that controls the numerical diffusion.
